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Abstract 

Two wavelength interferometry can in principle be used to measure changes in both 
temperature and concentration in a fluid, but measurement errors may be large if the fluid 
dispersion is small. This paper quantifies the effects of uncertainties in dn/dT and dn/dC 
on the measured temperature and concentration when using the simple expression dn = 
(dn/dT)dT + (dn/dC)dC. For the data analyzed here, ammonium chloride in water from - 
5 to 10 °C over a concentration range of 2 - 14% and for wavelengths 514.5 and 633 nm, 
it is shown that the gradients must be known to within 0.015% to produce a modest 10% 
uncertainty in the measured temperature and concentration. These results show that great 
care must be taken to ensure the accuracy of refractive index gradients when using two 
wavelength interferometry for the simultaneous measurement of temperature and 
concentration. 

Introduction 

Interferometry is an excellent method for measuring refractive index changes in fluids 
with high accuracy. When the refractive index changes are caused by changes in a single 
flow variable such as temperature, it is straightforward to convert interferometric data to 
engineering data. When two properties change the refractive index simultaneously, two- 
wavelength interferometry has been used to get an extra equation to solve for the 
additional unknown. 1 Similarly, three wavelengths can be used to solve for three 
unknowns. There are several sources of error in two wavelength interferometry. Some 
are common to interferometers in general, such as errors caused by fringe measurement 
error 2 ' 3 and errors caused by refraction and end effect. 4 However, as mentioned by Vest 5 
with regard to gases, two wavelength interferometry is very sensitive to small errors when 
the fluid is weakly dispersive. This paper presents an uncertainty analysis that quantifies 
the relationship between calculated temperature and concentration and the values used for 
dn/dT and dn/dC for the following relationship: 6 

An = (dn/dT)T + (dn/dc)c ( 1 ) 

where n is the refractive index, T is the temperature, and C is the concentration of the 
fluid mixture. Data specific to ammonium chloride in water is then analyzed because of 
its wide use for simulating superalloy solidification 7 in microgravity studies. It is shown 
that for ammonium chloride over the temperature range -5 to 10 °C over a concentration 
range of 2 - 14% and for wavelengths 5 14.5 and 633 nm the values of dn/dT and dn/dC 
must be known to high accuracy over the experimental range to ensure reasonable 
uncertainties in the final measurements of temperature and concentration. 
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Procedure 

Interferometry relates measured interference fringe spacings to refractive index 
differences as follows: 


A AM - 


= L- 


A n i 


( 2 ) 


A.v Av 1 

where A is the wavelength of light, L is the physical pathlength through which the light 
passes. As is a pathlength difference in the plane normal to the line-of-sight, and AM and 
An are the number of fringes, and refractive index change along As, respectively. The 
subscripts A refers to wavelength dependent quantities. The pathlength L is known from 
the experiment, the wavelength is known by choice of laser operating line, and AM/As is 
measured from an interferogram. Knowledge of these parameters permits the calculation 
of An/As for a particular wavelength. Rewriting equation ( 1 ) to describe the line-of-sight 
integration of interferometric measurements yields: 


dn | dn i dT dn i dC 
~ds'* = ~dT' x ~ds+~dC' k ~ds 


(3) 


In this case there are two unknowns, dT/ds and dC/ds. Therefore two wavelength 
interferometry should be sufficient for producing two equations required for a solution. 
Equations (2) and (3) are combined for each wavelength: 


A. dM, 
L ds 


dn | 

Hr 


dT dn _ , 
ds dC 


dC_ 
u ds 


(4) 


where i=l,2 represents the two wavelengths. Implicit in the assumption of a solution is 
that dn/dT and dn/dC are different for two different wavelengths, and that: 

(dn/dT ) | A1 (dn/dC ) | 


A1 


(dn/dT)\ u (dn/dC)\ u 

Forming a new variable from the experimentally determined quantities: 
A ,=(l,/LpM,/ds) 


(5) 


( 6 ) 


and shortening the nomenclature by the following new variables: 


dn | 

dn | 

Hr ' A ' 


dn | 

dn 

dC' A ' 



yields the simplified equations: 

A, = c lT (dT/ds) + t\ c [dC/ds) 

A 2 = c 2T [dT/ds) + c 2c ( dC/ds ) 

which can be used to solve for dT/ds and dC/ds: 
dT A A,c j£- 


(7) 


( 8 ) 


ds 

dC_ 

ds 


c )/ c :c C U C 2! 

A, c 2T - A 2 c ]T 

c ir c :i ~ c \t c 2c 


( 9 ) 
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Note that Equation (5) must be satisfied or neither of these quantities can be determined. 


What is the uncertainty in dT/ds and dC/ds (and therefore also in the measured 
temperature T and concentration C)? Assuming that uncertainties in the variables on the 
right hand side in Eqn. (9) are independent, uncertainty analysis yields: 8 



where e v represents the uncertainty in variable x. Uncertainties caused by the wavelength, 
geometric pathlength, and number of fringes per unit distance have not been included 
here. This analysis is meant to determine the effects of the refractive index variations 
with concentration and temperature only. From Eqn. (9) the partial derivatives can be 
calculated: 

d(dT/ds) ( -c 2C . W 

dc <^17^2 C' ^ I C^2T ) 

d(dT/ds) ( y u . \dT 

V C 1T^2C _ C iC C 2T / ^ 11 

d(dT/ds) c 2J R L /A, \ dT 

dc ^ c j ^ ^ic ^ i ) ds 

d(dT/ds ) c 2r RXJX, Y -iy W 

^ C |y ~ T ^ 2 C r ] R Ao j A] , ^ ^2C ) 

Similarly, 

djdCjds) ( c 2c - /?Ao/A, W 

dc jj’ ^ c j 7 '^* 2 c ^ 2 t C\jR Ao / A| ) ds 

d(dC/ds) ( c 2C Rh,/\ Y c l7 )dC 

dc ~>t \C \c^2T ^ 1 T ^2 C ^ 27 ^17^ ^2 / y ^ ^ 2 7' ^ (12) 

d(dC/ds) _ f -e 2r W 

* IC — C|j-C 2c - J (i.v 

d(dC/ds ) f c, r j <7C 

V ^ 1C^27 ^\T^2C J 
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where 


dMJds 

« S ^Ts’ ,l3) 

the ratio of the number of interference fringes per unit length measured from wavelength 
2 and 1 . From all of these equations, we can determine the effect of inaccuracy in dn/dT 
and dn/dC on the measured concentration and temperature. 

The variation of refractive index with temperature and concentration has been determined 
for two wavelengths 632. 8nm and 514.5nm by Spatz . 9 The reported relationship is: 

n HtNe = 1.3332 - (6.8276 ■ 1 <T 5 )t + ( 1.6764 • 10 _, )C 

n Ar = 1.3373 -(7.1 010- l(r 5 )7' + ( 1.7345- 10 _1 )c * * 4 ’ 

where temperature T is expressed in °C and concentration C is expressed in % by weight. 
These equations were stated to be valid for the a temperature range -5 °C to 10 °C and 
concentrations of 2% to 14% by weight NH 4 CI in water. Spatz also reports a value for 
R . 9 The required constants are summarized in Table 1. 


Symbol 

Value 


R = 

1.289 


XI = 

632.8 nm 


X2 = 

514.5 nm 

dn/dT @ XI 

C|T = 

6.8276E-5 (°C)‘' 

dn/dT @ X2 

C2T = 

7. 1010E-5 ( 0 cy‘ 

dn/dC @ XI 

Cic = 

1.6764E-3 

dn/dC @ X2 

C 2 C = 

1.7345E-3 


Table 1. Values of constants. 


Substituting these values into Equation (12) yields the following dependence of dT/ds and 
dC/ds on errors in the four constants c,j and c, c : 


' . = 396.8 — 

at /els ^ c, 


= -266.5 


A (dC/ds) 
dC/ds 


and 

Ac, r 

— ; — — — = zoj.j 

dT/ds ^ c, r 


= 396.8 


= 265.5 


= 266.5 


Ac if 
C 2T ) 

( A c 3c 

< C 2C 


A (dT/ds) 
dT/ds 


= —178.3 


= 178.2 


From Equations (15) and (16) it is clear that small errors in the four constants produce 
large errors in the calculated temperature and concentration distributions. In particular, a 
1% uncertainty in dn/dT for either wavelength produces a nearly 4-fold change in the 
calculated concentration, and a factor of 2.6 error in temperature. A 1% uncertainty in 
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dn/dC produces a factor of 2.6 error in the calculated concentration and 180% error in 
temperature. Note that using a value of dn/dT = 6.853E-5 instead of 6.8276E-5 for 
A,=632.8nm violates Eqn. (5) and prevents the determination of either concentration or 
temperature. This change in dn/dT represents a change of only -0.372%. Figure ( 1 ) 
shows n plotted as a function of temperature using both dn/dT = 6.853E-5 and 6.8276E-5 
(°C)' 1 in Equation ( 14a). The data are practically indistinguishable, certainly within the 
standard deviation in the data presented in Ref. 9. Yet this small difference produces 
huge changes in the measured concentration and temperature. 


These errors have been considered independently. Equation (10) determines to total 
uncertainty if each of the four constants can be considered to vary independently with 
equal weight. Normalizing Equation (10) yields: 


-i rrds 


cTT/ 


ds 


£jCd\ 


dC/ 



7 ds. 


and using the values from Table 1 yields: 


.,(17) 


(IT /ds 


dC/ds _ 


1 
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+ 

£ IT 

2665-^ 

+ 

£ ... 

- 1 78.3— L±L - 
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C \T _ 


C 2T . 


% . 



e, T 
398.6-^ 

+ 

£ .1 -r 

396.8 — ^ 

+ 

-2665 — 

2 

+ 

C\T . 


C 2T . 


% _ 



1.78.3 


265.5 


£ ij 2T 

c 2( 

£,;r 


(18) 


Assuming just a 0. 1 % error in each of the four constants yields: 
e m =0A529(dT/d S ) 

£ ICJx = 0.6756 (dC/ds) 


(19) 


that is, a 0. 1 % uncertainty in each refractive index variation produces a 45% uncertainty 
in temperature and a 67% uncertainty in concentration. These values scale linearly with 
assumed refractive index gradient uncertainties. To measure temperature and 
concentration to within 10%, dn/dT and dn/dC must be known to within 0.015%. 


Conclusions 

Uncertainty analysis was used to determine the accuracy with which dn/dT and dn/dC 
need to be made to ensure sufficient temperature and concentration accuracy from two- 
wavelength interferograms. For the data analyzed here, namely a temperature range from 
-5 to 10 °C and a concentration range of 2% to 14% by weight NH 4 C1 in water for 632.8 
and 514.5nm wavelengths, and using the simple linear expression relating refractive 
index to temperature and concentration described by Equation 1, the uncertainty must be 
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0. 015. or less for both dn/dT and dn/dC to measure the temperature and concentration to 
within 10%. This excludes all other error sources, including fringe measurement errors 
which can contribute up to 20% error. 2 Values for dn/dT are commonly obtained by 
using dn/dT = (-3/2)[n(n 2 -l)/(2n 2 +l)](3, where (3 is the coefficient of thermal expansion, 
but this expression is only accurate to about 2%. 10 Refractive index measurements 
accurate to four decimal places can be made by placing a liquid in a parallel sided glass 
cell, focusing a beam first on the front surface of the liquid and then the rear surface and 
measuring the focus shift. 1 1 Obtaining this accuracy requires distance measurements on 
the order of 100 nm and angular measurements on the order of a few arc seconds. 
Certainly the measurements must be made for the specific fluid under study over the 
specific temperature and concentration ranges, not extrapolated from published data. 

These very large final errors from small errors in dn/dT and dn/dC are the direct 
consequence of the slow variation of dn/dT and dn/dC with wavelength in ammonium 
chloride. This behavior is typical of fluids and gases. 

Implicit in these results is the assumption of the relationship between refractive index, 
temperature and concentration given by Equation 1. Recently researchers have refined 
this relationship to include quadratic terms 12 and fourth-power plus cross terms in the 
Cauchy equation. 7 This latter expression is complex, including thirty coefficients. 
Unfortunately, increasing the number of coefficients neither compensates for 
measurement noise nor increases the dispersion of the fluid and therefore should not 
reduce the sensitivity of two-wavelength interferometry to uncertainties in refractive 
index gradients. 
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FIGURE 1. Variation of refractive index with temperature for X = 632.8 nm, 
concentration ranging from 2% to 14% in 2% intervals, and two values of dn/dT: 
6.853E-5 and 6.8276E-5 (°C) '. 
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